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Abstract 

Proton beam transport in the context of proton driven Fast Ignition is usually assumed to be 
stable due to protons high inertia, but an analytical analysis of the process is still lacking. The 
stability of a charge and current neutralized proton beam passing through a plasma is therefore 
conducted here, for typical proton driven Fast Ignition parameters. In the cold regime, two fast 
growing Buneman-like modes are found, with an inverse growth-rate much smaller than the beam 
time-of-flight to the target core. The stability issue is thus not so obvious, and Kinetic effects 
are investigated. One unstable mode is found stabilized by the background plasma protons and 
electrons temperatures. The second mode is also damped, providing the proton beam thermal 
spread is larger than 10 keV. In Fusion conditions, the beam propagation should therefore be 
stable. 

PACS numbers: 52.57 Kk,52.59.-f,52.35.-g,52.35.Qz 
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Introduction - The Fast Ignition Scenario for Inertial Fusion has been proposed more 
than one decade ago as a solution to relax the symmetry constraints on the target irradiation 
in conventional schemes . Such relaxing allows for a lower laser energy, hence higher 
gain. The first Fast Ignition concept presented implied the propagation of a relativistic 
electron beam from the border of a pre-compressed target, to the dense core. It was quickly 
recognized that beam transport stability would be an issue for this scenario 3|, \^ , and many 
efforts have been devoted to document this point 

mm 

. Proton driven Fast Ignition was 
ater suggested as a probably more stable alternative, due to the higher inertia of the protons 
qI, . The scheme was made possible by the recently discovered possibility of generating 
intense beams of high energy protons, from the rear surface of solid targets irradiated by 



ultra-intense lasers 111 
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14| . This option has been largely investigated since then (see 



Ref. [lOj], and References therein), and alterative scenarios using more than one beam have 
already been presented [is], 16 1. 

Although a proton beam should definitely beam more stable than an electron beam, its 
stability within the context of proton driven Fast Ignition as not been addressed analytically 
so far. Further more, some simulation of the transport of a charge and current neutralized 
proton beam in a dense plasma, suggested that beam-plasma instabilities could also be an 
issue 17|. The goal of this paper is to clarify this problem by working out the kinetic 



theory of the proton beam instability within the alluded scenario. We thus start examining 
a cold model which shows that the e-folding time of the fastest growing mode is indeed 
much smaller than the beam time-of-fiight from the target border to its dense core. In this 
respect, unstable modes are driven by the interaction of the co- moving electrons with both 
the beam and the background protons. We thus implement the kinetic theory of the problem 
accounting for the temperature of the three species, and show that for realistic temperature 
parameters, the beam propagation should be stable. 

Cold Model - We consider the propagation of a proton beam of density and velocity 
Vb in a plasma. The plasma proton and electronic densities are Up and rie respectively. The 
beam is initially both current and charge neutralized, so that ui, + Up = Ue and UbVb = 
rigfe- Unlike the case of electron driven Fast Ignition, where protons are usually assume at 
rest due to their higher inertia, the background protons need here to be allowed to move, 
since the proton beam defines a slower dynamic. Given the typical beam energy of 15 
MeV considered in the scenario, protons are not relativistic. A non-relativistic formalism is 
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therefore implemented in the sequeL Such formahsm allows us to focus on the flow-aligned 
instabilities instead of having to search the full 2D unstable spectrum (see discussion in 
the Conclusion). The cold dispersion equation for flow-aligned perturbations varying like 
exp{ikz — iujt) reads 

l = ^+ + (1) 



with, 

r, 47rn„e^ , 47rn;,e^ n 47rnee^ 

where M and m are the proton and electron masses respectively. We now introduce the 
dimensionless variables, 

uj „ kvb m rib 

x = —, Z = , R=Y7^ a = — , (3) 

UJe UJe Jm Up 

SO that the dispersion equation ([1]) reads, 

R 1 Ra 1 1 

\ l+a I 

Figure [1] displays the growth-rate as a function of Z for a = 0.1 of the most unstable mode 
for a given Z . The local growth-rate extremum at low Z stems from the interaction of the 
co-moving electrons with the proton beam. The maximum growth-rate is found here for 
Z ~ 1 and X ~ 1 with, 

^ = ^(^«)'^'^- (5) 

With x/Z = uj/kvb ~ 1, this mode is therefore resonant with the beam. The all-spectrum 
fastest growing mode arises from the Buneman interaction of the electron beam with the 
background protons 13]. It is located at Z ~ 1/a and x ~ 0.05 with, 

^ = ^,R"'-^e. (6) 

With R = 1/1836, this Buneman mode is found growing at 5 = 0.052ci;e. For an initial 
electronic density rie = 10^^ cm~^, we thus have an e-folding time Tg = 1/6 = 10~^^ s. This 
e-folding time needs now to be compared with the time the proton beam takes to reach the 
target core. Assuming the beam needs to travel 100 fim at Vb ~ 0.146c (15 MeV), we find a 
transit time ~ 2 x lO^^^s, so that Tg/r^ > 200. This simple calculation of the growth-rate 
shows that the instability could reach 200 e-folding times before the beam reaches the target 




FIG. 1: growth-rate given by the cold dispersion relation ([3]) as a function of the reduced wave 
vector Z for a = 0.1, and R = 1/1836. 



core. We therefore turn to the kinetic theory in order to evaluate how kinetic effects slow 
down the instability. 

Kinetic Theory - In line with the Dawson's picture of the multi-beam model [2^ , Eq. ([1]) 
can be directly adapted replacing the terms corresponding to each species by their kinetic 
counterpart. For example, the proton beam of density Ub and velocity Vb is replaced by an 
infinite numbers of "micro-beams" of density nbfb{v)dv and velocity v, where fb{v) turns to 
be the beam velocity distribution function. The substitution operated is thus, 

1 r h{v)dv 



kvbf 



with 



h{v) 



M 



2nkBTb 



exp 



(tu — kvy ' 

2kBTb 



(7) 



(8) 



After substitution of the three terms of Eq. ([T]), the kinetic dispersion equation is found as. 



= 1 



p W 



+ 



2u^ 
- W 



k'^v^ 



Te 



VTp 

uj/k — V, 

VTe 



k'^v^ 



Tb 



uj/k — Vb 



VTb 



where. 



te-* dt 



and the three thermal velocities read, 

2 2kBTp 



-'Tp 



M 



2kBTb 
M 



2kBT, 



m 



(9) 



(10) 



(11) 
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Equation ([9]) can be of course equally derived from the linearized Maxwell- Vlasov system of 
equations. In terms of the dimensionless variables defined in Eqs. ([3]), the equation reads, 



= 1 + 



+ 



R 



1 



Ra 



rp- 
1 



W 
W 



X 



ZpTp 

x/Z- 



PTb 



+ 



1 



with, 



PTp 



VTp 



w 



x{l + a)/Z — a 
PTe(l + a) 



(12) 



pTb 



(13) 



VTb VTe 
, PTe = • 

Vb Vb Vb 

At this junction, we consider typical proton driven Fast Ignition parameters, namely a quasi 
mono-energetic proton beam of mean energy 15 MeV, energy spread 1 MeV and density 
10^° cm~^, entering a plasma of density 10^^ cm~'^ and temperature 1 keV for both protons 
and electrons. These values set the dimensionless parameters defined by Eqs. (13|T3|) to 



a 



0.1, = 8.10-3, PTb = 0.25, PTe = 0.35. 



(14) 



The three distribution functions are displayed on Figure O allowing for a intuitive evalu- 
ation of kinetic effects upon the two aforementioned unstable modes. The large k mode, by 
virtue of its small phase velocity, is expected to interact mostly with the electrons and the 
background protons thermal spreads. For this mode, the beam is virtually cold. 

As evidenced on the plot, the small k mode is beam resonant, but the large electronic 
thermal spread should significantly affect it. For this mode, the background protons are 
almost cold, but the thermal spreads of the remaining species should be important. 

These trends are now confirmed by Figures [3] and HI where the numerical evaluation of 
the growth-rate has been performed for both modes in terms of the reduced wave vector 
Z = kvfj/ujp, and for various temperature parameters. From the numerical point of view, 
the full kinetic dispersion equation (fT2l) has been solved in each case, using the "FindRoot" 
routine of Mathematica and providing the cold solution as an initial guess for the algorithm. 
Figure [3l documenting the small k growth, has been plotted setting the background protons 
temperature to its realistic value pTp = 8.10"^ in every case (except the cold one). As 
expected, this mode is thus found virtually insensitive to this thermal effect. But both 
the beam and the electron thermal spreads contribute here to stabilize the waves. For a 
beam temperature parameter pxb = 0.25, the mode is almost stabilized with pTe = 0.1, 
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FIG. 2: (Color online) Distribution functions of the three species involved. The positions of the 
phase velocities of the two most unstable modes are located. 



FIG. 3: (Color online) growth-rate of the low k unstable mode as a function of Z = kvb/ujp, and 
for various temperature parameters. 

which remains much smaller than the expected value. Complete stabilization is achieved 
for pTe = 0.24 (Tg = 0.47 keV), so that the part of the spectrum is completely stabilized 
by PTe = 0.35 (Te = 1 keV). The stabilization of this mode is thus provided by the beam 
and the electron thermal spread. Setting every temperature parameters to their realistic 
value, and varying the beam one, it is numerically found that stabilization is achieved from 
PTb > 2.6 X 10^^. This means that a beam with a thermal spread smaller than ~ 10 keV, 
far smaller than the 1 MeV observed so far, would not propagate stably. 

The evolution of the large k unstable mode is pictured on Figure HJ where the beam tem- 
perature parameter is set to its realistic value in any cases, evidencing the small influence 
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FIG. 4: (Color online) growth-rate of the large k unstable mode as a function of Z = kvb/ojp, and 
for various temperature parameters. 

it has on overall thermal effects. For these short wavelength Buneman modes, stabiliza- 
tion comes from the background protons and the co-moving electrons thermal spread. For 
background protons realistic temperature parameters, stabilization is almost achieved at 
PTe = 3.10^^. Complete stabilization is numerically found from pxe > 3.3 x 10~^ (Tg = 10 
eV), which is far lower than the expected thermal spread. It is important to note that the 
stabilization of this mode does not rely on the beam quality, and is only provided by the 
initial plasma temperature. 

Conclusion - The stability of an intense proton beam passing through a plasma has 
been analyzed. In line with previous numerical approaches pJZj], the beam is assumed to 
be both charge and current neutralized. A simple cold model describing the flow aligned 
unstable modes evidences two Buneman-like unstable modes, arising from the co-moving- 
electrons/proton-beam, and CO- moving-electrons/background-protons, interactions. Note- 
worthily, the fastest growing mode in the cold regime is found to have an e-folding time 2 
orders of magnitude lower than the beam time-of-flight to the dense core. A kinetic calcula- 
tion is then undertaken, which shows that for realistic temperature parameters of the three 
species, flow-aligned instabilities should deflnitely be stabilized. 

As previously mentioned, we focused here on flow-aligned unstable perturbations. When 
it comes to the unstable modes involved in the case of a relativistic electron beam for ex- 



amp 
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e, it is now known that Filamentation or Oblique modes can govern the linear phase 
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23|, with a wave- vector not parallel to the beam direction. Such shifting of the 



dominant mode is mainly due to anisotropic relativistic effects which make the particles 
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"heavier" to move along the beam axis than in the perpendicular direction. In the non- 
relativistic regime examined here, the unstable spectrum is still 2D, but the growth of 
non-parallel instabilities is hindered by thermal effects which prevent (or simply slow-down) 
the formation of filamentary structures Q]. Electrostatic modes growing along the beam 
direction are thus expected to be the most unstable ones. It could happen that once these 
modes have been stabilized, the system remains unstable with respect to the Filamentation 
instability. Indeed, situations have been evidenced where a relative drift still excites the 
later, while the formers are stabilized (2^. At any rate, the growth of Filamentation must 
be slower in the kinetic regime than in the cold one. The e-folding number of the Filamen- 
tation instability for the present case has thus been checked, and found lower than 3. This 
instability should thus not be a threat in the present case. 

Finally, we here considered an infinite and homogenous plasmas, while the realistic sit- 
uation involves a density gradient. In this respect, the uniform plasma approximation is 
valid providing the wavelengths implied in the instability are short compared to the scale 
length of the density gradient. This question has been addressed in a previous work {25]. On 
the one hand, the wavelengths relevant to the unstable spectrum are typically the plasma 
skin depth which is 0.17 fim aX n = 10^^ cm^'^. On the other hand, an exponential density 
gradient from the target border to its core sets a gradient scale length A ~ 10 fim. The 
present approximation is thus valid from the onset of the beam transport, and even more 
justified at later time since the plasma skin depth is monotonically decreasing as the plasma 
density rises. 
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